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Abstract. The current popular cosmological models have granted to the population of dwarf satellite galaxies 
' a key role: their number, location and masses constrain both the distribution of dark matter and the physical 

evolution of their hosts. In the past years, there has been increasing observational evidence that objects with 
masses of dwarf galaxies can form in the tidal tails of colliding galaxies and speculations that they could become 
' satellite-like galaxies orbiting around their progenitors and thus be cosmologically important. Yet, whether the 

so-called "Tidal Dwarf Galaxies" (TDGs) candidates are really long-lived objects and not transient features only 
present in young interacting systems is still largely an open question to which numerical simulations may give 
precise answers. We present here a set of 96 N-body simulations of colliding galaxies with various mass ratios 
and encounter geometries, including gas dynamics and star formation. We study the formation and long-term 
evolution of their TDG candidates. Among the 593 substructures initially identified in tidal tails, about 75% 
fall back onto their progenitor or are disrupted in a few 10^ years. The remaining 25% become long-lived bound 
(^-^ . objects that typically survive more than 2 Gyr with masses above 10^ Mq. These long-lived, satellite-like objects, 

' are found to form in massive gaseous accumulations originally located in the outer most regions of the tidal tails. 

\^ , Studying the statistical properties of the simulated TDGs, we infer several basic properties that dwarf galaxies 

■ should meet to have a possible tidal origin and apply these criteria to the Local Group dwarfs. We further found 

that the presence of TDGs would foster the anisotropy observed in the distribution of classical satellite galaxies 
^ ^, around their host. Identifying the conditions fulfilled by interacting systems that were able to form long-lived tidal 

I ' dwarfs - a spiral merging with a galaxy between 1/4 and 8 times its mass, on a prograde orbit, with an orbital 

O I plane inclined up to 40 degrees to the disk plane - and estimating their fraction as a function of redshift, we 

roughly estimate their contribution to the overall population of dwarfs. We conclude that a small but significant 
^ ' fraction of them - a few percent - could be of tidal origin. This number may be underestimated in particular 

environments such as the vicinity of early type galaxies or in groups. 
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1. Introduction studies that dwarf galaxies have a cosmological primor- 

■ r 1 r 11- 1 • 1 dial origin and were not formed recently. However, there 

The properties or dwarf satellite galaxies that surround . . . , _li _l i ■ ^ -^1 ^ ■ i r 

f ^ , is growing evidence that objects with masses typical or 

spirals and ellipticals are oiten used to constrain cosmolog- , r i ■ r • i-j i ^ -i ^ j i 

. ^ . Ill • 1 J 1., J • dwari galaxies form m the tidal tails th at surround col 



ical models. Their number has been actively debated since ^j- j ■ ■ i i • / ^^^ paT T~ 

, , , , , liding and merging spiral galaxies (e.g. l Uuc fc Mirabe 
cosmological simulations have shown that the low-mass ' ^ ' ' — — ■ ' — 



dark halos are much more numerous around their mas- 
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T993lHibba,rd et a,lJl2nmUMendes de Oliveira, et ~ 
Temporin et alJbOOa iKnierman et al.ll2003l) . Many 
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, „^ , r— — — "—11 I r— -I I , I 01 these so-called iidal Dwari Ga laxies (IDCjs) ap 

Milky Way (e.g. |Klypm et al.|[l99i |Moore et ajJllMi. , ..if_„...u..i.„ .Ki..f. \ nJJ:,.A a 



, p „. pear to be self-grav itating objects llBraine et al.1 12OO0I: 

Moreover, the spatial distribution 01 satellites can pro- fS j 7 1 lonn/il\ Tfii, i-j 1 u- j. 1 r j 

_ . „ ' . ^ , , , r 1 1 1 Y lljournaud et alJ2004l) . It these tidal obiects are long-lived. 

vide mtormation on the mass and shape of dark haloe s ^ -u ^ ^ j- i 1 w r j c 

— : — '■, n k 1 — : irrrm rr^ — 1| „„„ they could contribute to the total population ot dwart 

llZaritskv et alJll997t lBrainerdll200a iLee k Kanell2006l) . , ^ ^ r rr.. ■ . 

_r ~ — ~ ' — 7"^ ~" — ~^ ~^ 7 satellites m addition to primordial dwarts. iheir statis- 

Fmally, they play, as building blocks, a maior role on the , . , , . , ,-n j 1.1 -i 1 

. ' , . r . . . , , . , tical properties would then be modihed and the possible 

physical evolution of massive galaxies with which they .-^ 1-1 11 iii ,.1 

^ „ T • n 1-1 constraints on cosmological models would have to be up- 

will eventually merge. It is generally assumed 111 such (^jg^^^^jj 
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may cause a rapid orbital decay. iHibbard fc Mihoj l)l995|) 
have shown that a large part of the once expelled tidal 
material falls back onto the parent spiral galaxies in a few 
hundreds of Myr, but the time-scale increases to Gyr in 
the outer regions; thus the tidal dwarfs formed there may 
have a significant lifetime. Also, dwarf galaxies may be 
disrupted by the tidal field of their progenitor even if the 
internal orbits of t heir stars are not res onant with their 
orbital period (e.g. iFleck fc KuhiJEoO^ : this disruption 
process can take a few billion years depending on their 
mass, orbit, and concentration. Thus, the lifetime of tidal 
dwarfs is a priori rather uncertain and their cosmological 
importance far from being proven. 

Actually, no real old Tidal Dwarf galaxies, still sur- 
viving after the merger of their progenitors and the van- 
ishing of the umbilical cord linking them to their par- 
ents, has yet been unambiguously fou nd although candi- 
dates were discus sed in the literature l|Hunter et al ] l20nflt 
IDuc et a'l]l2004bL and references therein). Because of the 
difficulties to identify a tidal origin in evolved galaxies, 
numerical simulations appear to be unavoidable to pre- 
dict their number. Sofar, numerical modeling has mainly 
been used to study the formation of objects in tidal tails. 
The gravitati onal collapse of p re-existing cl ouds or regions 
of tid al tails l Elmegreen et al.t. 199^1 .Barnes fc HernauistI 
can lead to the formation of clumps along tidal 
tails with typica l m asses of 10^~^ Mp. A ccording to 
IWetzsteinI l)2005l) and IWetzstein et all 1)2006^ . a massive 
gaseous component is required to achieve this collapse. 
More massive accumulations of matter can form in the 
outer regions of tidal tails when dark halos are assumed 
to extend much further than the optical radius of galaxies 
ijBournaud et al.ll2003j) . Their formation is initiated by a 
kinematical mechanism, and self-gravity can make these 
accumul ations of matter collapse into dense star-forming 
objects tone et al.ll2004al hereafter D04). Simulations as- 
suming massive and extended gaseous disks in the collid- 
ing spirals pr oduce objects th at seem to become satellites 
galaxies (e.g. lCox et al.ll2'004i D04) but the conditions for 
their survival have not yet been studied in detail. The fate 
of classical infalling satellites has been studied i n numer- 
ous simulations fe.g. iMaver et al.l l2001a"b'. 2002). In such 
models, the dwarfs were assumed to contain a somehow 
"shielding" layer of d ark matter, while tidal dwa r fs do not 
in a C DM scenario (|Barnes fc Herna uistI Il993 ). iKroupal 
studied the survival of tidal-like dark-matter free 
dwarf satellites but started his simulations once the galax- 
ies had already been formed and hence could not address 
the role of their initial distribution on their evolution. He 
argued that the objects produced in his model - with final 
masses of only 10^ Mq - become long-lived after having 
lost up to 99% of their initial mass. These results cannot 
be easily extended to more massive objects since the latter 
undergo more dynamical friction and may have a differ- 
ent morphology. Hence, self-consistent simulations of the 
formation and subsequent evolution/survival of TDGs are 
still required to estimate the contribution of such objects 
to the population of dwarf satellite galaxies. 



In this paper, we study the long-term survival of mas- 
sive dwarf-like objects formed in tidal tails. A set of 96 
numerical simulations of galactic encounters were carried 
out, varying the geometrical parameters of the collision 
and the relative masses of the colliding galaxies. They are 
used to predict the number and distribution of long-lived 
- with life expectancies greater than one billion year - 
tidal objects that will eventually look-like classical dwarf 
satellites. The numerical techniques and parameters are 
described in Sect 2. In Sect. 3, we detail our results on the 
formation criteria and survival time of the tidal objects; 
we then determine their characteristics once they have 
become long-lived satellite galaxies. In Sect. 4, we com- 
pare them with the properties of respectively the currently 
forming TDG candidates observed in young interacting 
systems and the satellites around massive hosts such as 
those identified in the Sloan Digital Sky Survey and in 
the Local Group. In particular we discuss rough calcu- 
lations on the expected fraction of dwarfs of tidal origin 
among all classical ones. Our conclusions are summarized 
in Sect. 5. 



2. Numerical simulations 

2.1. Numerical techniques and parameters 

The simulations presented in this paper use the same 
particle-mesh FFT code as in D04. In particuler, a sticky- 
particles scheme is used to model the ISM^ and star for- 
mation is described by a local Schmidt Law. Initial condi- 
tions for disk galaxies and dark haloes are largely similar 
to those described in D04. Galaxies are assumed to con- 
tain 15% of gas (gas-to-visible mass ratio, including gas 
outside the stellar disk). Each disk galaxy is embedded 
in an extended halo that maintains a flat rotation curve 
up to ten times the stellar disk radius, consistent with 
cosmological predictions. 

We use a 512'^ cartesian grid to compute the gravita- 
tional potential, with a gravitational softening of 380 pc. 
The most massive galaxy is modeled by 10® particles for 
each component (stars, gas, dark matter). Its stellar mass 
is 2 X 10^^ Mq and it contains 10% of gas. Its stellar 
disk radius is 15 kpc and its gaseous disk radius 45 kpc 
- HI is indeed observ ed much beyond the opt ical radius 
in spiral galaxies (e.g. lRoberts fc Havneslll994j) . The dark 
halo is modeled by a pseudo-isothermal sphere, with an 
asymptotic velocity of 205 km s^^ and a core radius of 
10 kpc. The circular velocity at the dark halo virial ra- 
dius (280 kpc) is 220 km s~^. All the timescales discussed 
in this paper are hence realistic for galaxies similar to 
the Milky- Way and must be rescaled accordingly for other 
galactic mass/virial velocities. The number of particles in 
the companion depends on its mass; its size is scaled by 



^ D04 explained that the modeling of the ISM is not crit- 
ical for the form ation of massive tidal objects, and indeed 



ICox et al] ll2004f) have obtained the formation of comparable 
objects with an SPH code 
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the square root of its mass to keep the stellar disk surface 
density constant. 

Ninety-six galaxy mergers with various mass ratios and 
various orbital parameters have been simulated. The de- 
tailed parameters for each run are given in Appendix 1X1 
The range of values explored for each parameter is sum- 
marized in Table In order to increase the statistical 
size of our sample of TDG candidates, we have privileged 
cases favorable to the formation of long-lived tidal dwarfs. 
For instance, because retrograde orbits are found not to 
form long-lived TDGs, we have mainly simulated prograde 
cases. Of course, we take into account this bias in our sta- 
tistical study. 

The simulations have been run until 2 Gyr after the 
first pericenter of the relative orbit of the two galaxies. 
The " merger" typically occurs 50-300 Myr after this peri- 
center. We were actually able to follow the evolution of 
the objects formed in tidal tails for more than 1.5 Gyr 
after their formation which is enough to know whether 
they have become long-lived dwarf satellites, i.e. objects 
orbiting around their massive host for a few dynamical 
times. 

2.2. Detection of massive substructures in tidal tails 

We analyzed the spatial distribution of matter at different 
epochs of the simulations. We identified the massive sub- 
structures along the tidal tails that may become long-lived 
"Tidal Dwarf Galaxies". The tidal objects were sought in 
the gas component plus the young stars formed in situ af- 
ter the periaster. Old stars from the parent's disks were 
not considered as they seem to play a minor role in the 
formation of structures in tidal tails, at least in the most 
recent scenarii put forward: the kinematical pilling-up of 
gas at the tip of a tail (D04) or the dynamical forma- 
tion of gravitational cl umps in massive gaseous tidal tails 
(Wetzsteiu et h]M()^. 

The criterion for an "object" to be detected is that 
the mass of gas and young stars within a diameter d must 
be larger than a threshold Mq. Two combinations of mass 
threshold / maximal diameter have been used : 

Mo = 3 X 10** within d < 6 kpc 

or 

Mo = 1 X 10® within d < 3 kpc 

We do not use the combination 

Mo = 1 X 10** within d < 6 kpc 

since it was sometimes found to include parts of tidal 
tails that visually do not correspond to any substructure. 
Objects of 10® Mq contain 2 x 10^ to 10^ particles, de- 
pending of their gas fraction, so that objects above our 
mass threshold contain enough particles to be well re- 
solved. 

The identification of objects was made around pixels 
with a density larger than /io = Mg/{Tr(P /4), and the 
center of the spherical boundary (the 6 or 3 kpc diameter) 



was computed to maximize the included mass. At this 
stage, one object may in fact correspond to two different 
condensations of matter. We decided that one object must 
be separated in two if two components more massive than 
the mass threshold and separated by a density fi < fio 
were found. Note that we have rarely observed objects that 
were first considered as double to merge into a single one, 
so that the somehow arbitrary threshold used to separate 
them was realistic. 

The mass center of an object, used to determine its 
" position" , was defined as the mass center of pixels with a 
density larger than /io within the object external radius. 
We only considered objects located outside 1.5 times the 
initial stellar disk radius of each parent spiral galaxy, in or- 
der not to include spiral arms or other features belonging 
to the outer disk. Objects located at smaller radii would 
anyway not survive for a few 10® yr before merging and do 
not deserve being considered as potential TDGs. In all the 
paper, t = corresponds to the first pericenter between 
the two parent spiral galaxies, not to the beginning of the 
simulations. The criterion for the detection of massive sub- 
structures in tidal tails was first applied at t = 200 Myr. 
Later on, every 100 Myr up to f = 2 Gyr, each simulation 
was re-analyzed to know whether detected objects were 
still present (i.e., they were still obeying our mass/size 
and radial distance criterion) and identify new ones that 
were not formed at earlier times. 

An example of the detection of sub-structures in a sim- 
ulated merger, using our criteria, is shown in Fig. ^ 

Our detection algorithm based on a density threshold 
might be considered rather basic compared to more sophis- 
ticated algorithms, like the friend-to-friend one. Another 
approach could also have been to select only the grav- 
itationally bound structures that were likely to become 
long-lived dynamically independent objects. However, the 
goal of our simulations is mainly to make relevant com- 
parison with observations. Projection effects and a limited 
resolution make the identifica tion of bound sub-struc tures 
in real systems difficult ( see iBournaud et al.l l|2004j) and 
iHibbard fc BarnesI (|200J). Most of the so-called TDG 
candidates described in the literature were actually identi- 
fied based on their shape rather than on kinematical data. 
Because we wished to treat our simulations as would be 
real images, we decided to use a criterion solely based 
on the apparent morphology. We thus first included some 
unbound objects, like observations do, but following their 
evolution, we could a posteriori reject them. Beside, this 
allowed us to study the process of tidal disruption and the 
falling back of tidal material. 

3. Results: statistics on the simulated merging 
systems 

3.1. Number, lifetime and masses of tidal 
sub-structures 

Over our sample of 96 simulations (184 parent spiral galax- 
ies) , we identified 593 massive substructures in tidal tails, 
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Parameter 


simulated range 


range most favorable to long-lived TDG formation 


Relative velocity 


50 to 320 km s"^ 


50 to 250 km s"' 


Impact parameter 


15 to 200 kpc 


30 to 200 kpc 


Orbit inclination 


to 60 degrees 


to 40 degrees 


Orbit orientation 


Prograde/Retrograde 


Prograde only 


Mass ratio 


10:1 to 1:10 


4:1 to 1:8 



Table 1. Summary of the explored ranges for orbital and mass parameters, and indication of the range favorable to 
the formation of long-lived TDG. The detailed parameters, and number of long-lived TDGs formed in each simulation, 
are given in Appendix^ The relative velocity of the colliding galaxies is computed for an infinite distance, neglecting 
dynamical friction before the beginning of the simulation. The orbit inclination for each galaxy is the angle between 
the disk plane and the orbital plane. The mass ratio is given as the TDG-progenitor galaxy to tidally disturbing galaxy 
mass ratio. 



i.e. an average of 3.2 tidal substructure per parent galaxy. 
423 of these structures were first detected at i = 200, and 
most the other ones at t — 300 or t — 400. 

At t = 500 Myr, 207 objects - about one third of 
all the detected structures - still exist. Conversely, the 
two other thirds have disappeared in less than 300 Myr. 
At < = 1 Gyr, 143 objects are still present, and 119 at 
t = 2 Gyr (20% of the total). 63% of the vanishing objects 
have fallen back onto their parent galaxies; the remaining 
objects "disappeared" because their mass dropped below 
our 10^ Mq threshold. This can be due to tidal disruption 
or to the fact that they were not gravitationally bound. 

The distribution of their radius of formation (i.e. the 
radius at which they are first detected) along the tidal 
tails is shown on Fig. [21 It is roughly bimodal, with two 
thirds of the objects formed within the first half of the ra- 
dial extent of the tail, while the other third is formed near 
the tip of the tails. A close analysis of the simulations 
(e.g.. Fig. nj shows that the first kind of objects corre- 
spond to loose unbound structures and to gravitational 
clumps fo rmed all along the tails (su ch as the ones de- 



scribed bv ' Barnes & HernauistI 1 1992f ) or lWetzstein et al\ 
(j2006J), while the second kind of objects correspond to 
massive accumulations of matter formed near the tip of 
tails according to the mechanism described by D04. 

Each type of tidal objects appear to have different life- 
times: 

— self-gravitating accumulations of material formed near 
the tip of tidal tails generally survive for several 10^ yr 
or a few Gyr 

— other objects formed at smaller radii tend to disappear 
in less than 1 Gyr, and even often in only 500 Myr. 

Both categories of objects also show differences in their 
respective masses, as shown in Fig.O Most structures de- 
tected a.t t = 200 have a mass smaller than 10^ Mq; but 
those surviving at least t = 1 Gyr have the highest masses, 
with 45% of them exceeding 10^ solar masses. 

Therefore, according to our simulations, only the ob- 
jects formed at the tip of tidal tails can become long-lived, 
massive, dwarf satellite galaxies. The less massive objects 
that condensed at smaller radii along the tails tend to fall 
back rapidly on their parent galaxies or rapidly lose most 
of their mass, probably contributing to the population of 



super stellar clusters or even globular clusters observed 
around merger remnants. Our simulations do not have the 
resolution to study them in details. 
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Fig. 2. Number of objects as a function of their formation 
radius, defined as the radius at which they are first de- 
tected, normalized by the total radial extent of the tidal 
tail - the tidal tail extent is defined to contain 98% of 
the mass that has been pulled out 1.5 times the initial 
gaseous disk radius. The results are shown for the objects 
existing a.t t — 200 (green dotted), 500 (red dashed), and 
2000 Myr (blue solid). The longest-lived objects are those 
formed near the tip of tidal tails, many of which still exist 
after 2 Gyr. On the other hand, gravitational clumps and 
other tidal debris formed at smaller radii rarely survive 
more than 500 Myr. 



3.2. Constraints for the formation of long-lived tidal 
objects 

We now focus on the "long-lived" objects able to survive at 
least 1 Gyr - those preferentially formed in the outermost 
regions of interacting systems. The table in Appendix fXl 
indicates whether such objects are formed for each indi- 
vidual run and each parent galaxy. Table ^ summarizes 
the constraints obtained for the individual parameters of 
the collision (eventhough they are not fully independent). 
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Fig. 3. Mass spectrum of the objects detected in tidal tails 
at t = 200 Myr (green dashed) and t = 1 Gyr (red solid). 
The most massive objects that have the longest life are 
formed near the tip of tidal tails. 



The three main restrictions can be summarized as follow: 
a given spiral galaxy can form long-lived tidal objects pro- 
vided that: 

— the orbit is prograde: the angle 9 between the disturbed 
galactic dis k spin and t he com panion orbital spin (see 
Figure 8 in iDuc et alJ ll200(ll) 'l. is smaller than 90 de- 
grees. 

— its orbital plane is inclined by less than 40 degrees 

— the mass of the companion is at least one fourth of that 
of the target galaxy (otherwise the tidal field cannot 
drive matter far enough from the parent galaxy) and 
less than 8 times the target galaxy mass (otherwise 
tidal objects may form but rapidly fall onto the very 
massive companion). In other words, for mass ratios in 
the 1:1-8:1 range, at least one of the two spiral galaxies 
can form long-lived objects. For mass ratios in the 1:1- 
4:1 range, long-lived objects can form in the material 
of both galaxies. 

When these three conditions are achieved, most mergers 
will lead to the formation of long-lived tidal substructures, 
except those with small impact parameters (typically be- 
low 30 kpc), but the latter condition is a minor constraint 
since the small impact parameters are statistically far less 
frequent. 

When a spiral galaxy merges with a companion inside 
the parameter range described above, according to the 
results given in Table 1X1 on the average 1.9 objects are 
still visible at t = 500, and 1.3 survive at least 1 Gyr. 

3.3. Statistical properties of dwarf galaxy satellites of 
tidal origin 

3.3.1. Spatial 3D distribution 

As shown before, the long-lived tidal objects will behave 
as satellite galaxies, orbiting around their parent galaxies. 
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Fig. 4. Distribution of the long-lived TDGs (surviving at 
least 1 Gyr) in the radius-inclination plane. The radius 
is measured from the progenitor galaxy center; the incli- 
nation is the latitude in the progenitor disk frame. Data 
are accumulated from 500 to 2000 Myr. The orbital tracks 
of two objects are shown. They have the same initial ra- 
dius, but one (shown in red dashed) was formed on a more 
eccentric orbit attested by a larger radial excursion; it un- 
dergoes a larger radial decay. 
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Fig. 5. Latitudinal distribution of long-lived TDGs 
around their progenitor galaxies (red), compared to an 
isotropic distribution (green dashed). TDGs are prefer- 
entially distributed towards the equatorial plane of their 
progenitor disk. This distribution has been computed with 
data accumulated from t = 500 to 2000 Myr, for TDGs 
that survive at least 1 Gyr. The anisotropy of the distri- 
bution increases with radius, as detailed on Fig|S| 



We analyze in the following their spatial distribution, as 
derived from our simulations. To increase the size of the 
statistical sample, we have cumulated the data on the po- 
sitions from t ~ 500 to 2000 Myr for objects that survive 
at least 1 Gyr. The number of "detections" is then 1587, 
corresponding to 143 individual objects tracked at 6 to 
16 epochs of their evolution depending of their lifetime. 
Cumulating data not only enables to increase the statis- 
tical size of the sample, but also to make the comparison 
with samples of real colliding galaxies that are observed 
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at different times of their evolution. Of course, in the case 
of the simulations, not all the data points are statistically 
fully independent. 

We measure at each time the radius and inclination 
with respect to the parent galaxy disk. This inclination 
i is defined as the angle between (i) the parent galaxy 
disk, or, after the merger, the main flattening plane of 
the parent galaxy material^, and (ii) the line linking the 
parent galaxy mass center to the tidal dwarf mass cen- 
ter (Fig. d). We show on Fig. 01 their distribution in the 
(radius, inclination) plane. The inclination distribution is 
shown on Fig. El for all the objects, and on Fig. El for the 
objects separated into four quartiles depending on their 
initial radius (0-82, 82-121, 121-157 and >157 kpc). 

The distribution of the long-lived tidal objects around 
their parent galaxies (or the associated merger remnants) 
is strongly anisotropic. The fraction of objects found 
at galactic latitudes larger than 60 degrees is 4 times 
smaller than what is expected for an isotropic distribu- 
tion (Fig. jS}. Tidal dwarfs are mostly distributed in the 
equatorial plane of their progenitor disk. As mentioned 
earlier, they are preferentially formed when the orbital 
plane and the parent spiral galaxy plane are inclined by 
less than 40 degrees. In such cases, the tidal tails and the 
sub-structures formed within them are close to the disk 
plane. In a few rare cases, interactions with the progeni- 
tors and, potentially, with other tidal objects, drive some 
objects far from the parent disk plane. 

The anisotropy of the distribution increases with ra- 
dius, as seen on Fig. Indeed, tidal tails are more ex- 
tended in the close-to-coplanar cases and the objects form- 
ing at their tip are located at larger distances (Fig. ISJ. 
Finally, the most prominent tidal substructures have a 
more anisotropic distribution. We measure an average 
value of the inclination < i >= 31 degrees for objects 
less massive than 10^ Mq, and < i >~ 22 for the more 
massive ones. This, again, is explained by the fact that 
the most massive objects are formed on low orbital incli- 
nations (i.e. nearly-coplanar galactic encounters). 

3.3.2. Orbit eccentricities 

We determined the eccentricity of the orbits of the tidal 
satellites and followed its variation as a function of time. 
The "eccentricity" of the orbit of a given object at a given 
epoch was defined as follow: the orbit that this object 
would have if the potential of the merging system did 
not evolve is computed. A maximal and a minimal radii, 
respectively and r_ measured from the mass center 
of the merging/merged spirals, are found. The eccentric- 
ity e is then derived assuming that the orbit is in first 

^ In major mergers forming elliptical galaxies, the main flat- 
tening plane of the remnant is often found to correspond to the 
orbital plane, more than to the parent galaxy plane. However, 
since TDG-forming cases correspond to low orbital inclina- 
tions, the parent galaxy plane and the orbital plane are never 
very different. 
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Fig. 8. Radial distribution of the long-lived objects, de- 
pending on the inclination of the companion orbital plane 
with respect to the progenitor disk plane: orbital inclina- 
tions smaller than 27 degrees (blue solid) and larger than 
27 degrees (green dashed). Low-inclination orbits form the 
longest tidal tails and statistically produce TDGs at larger 
radii than more inclined orbits. 



order elliptical: e — (r+ — r_)/(r_|_ -|- r_). We show in 
Fig-Elthe distribution of e for the long-lived objects when 
they are formed, i.e. at their first detection (generally at 
t = 200 Myr). Nearly 90% of them were born on orbits 
with an axis ratio larger than 0.6; very eccentric orbits are 
rare. Objects with high eccentricities have a smaller peri- 
center radius, and are then exposed to a faster and more 
disruptive tidal field from the central galaxy; they also un- 
dergo a higher dynamical friction, causing them to quickly 
fall back onto their progenitor. At 2 Gyr, they have dis- 
appeared. The short-lived tidal structures, also presented 
in Fig. El had, without surprise, high eccentricities. 

We illustrate this point in Fig. ^ where we show the 
orbital track of two objects, initially formed at similar 
radii, in similar colliding systems (2:1 mergers, run 5 and 
6), but with different orbital eccentricities. The most ec- 
centric tidal dwarf is found to undergo a stronger orbital 
decay, under the effects of dynamical friction. It also loses 
18% of its mass from t = 0.5 to 2 Gyr, while the less ec- 
centric dwarf loses only 8% of its mass during the same 
period. 

3.3.3. Relative velocities 

Since they move on low eccentricity orbits, the tidal satel- 
lites generally have a velocity close to the circular veloc- 
ity of the parent gravitational system (dominated by the 
merged dark halo) . The resulting distribution of the veloc- 
ities at i = 1 Gyr is shown on Fig. ^1 It ranges between 
50 and 400 km s^^ (for a massive progenitor with a virial 
velocity of 220 km s^^). Assuming that systems are ob- 
served under random orientation, we computed the distri- 
bution of the projected linc-of-sight velocities (see Fie:. ll0() . 
The average absolute value is < IFlosI >= 115 km s^^ 
with respect to the mass center of the central merger rem- 
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Fig. 9. Statistical distribution of the initial eccentric- 
ity (at t = 200 Myr) of the objects able to survive 
at least 1 Gyr which became satellites (red dash-dotted 
histogram), and final eccentricity for those surviving at 
least 2 Gyr (green histogram) . The dashed blue histogram 
shows the initial eccentricity of objects surviving less than 
1 Gyr, renormalized to the same number of objects than 
for the red histogram. 



nant. In the case of young systems, a few 10^ yr after 
the merger, tidal tails a re more easily seen edge-on (e.g. 
iHibbard fc Barnesll2004l) so that the real line-of-sight dis- 
tribution is not isotropic. This can statistically decrease 
the projection angle, which in turns increases the observed 
line-of-sight velocities. The observed distribution around 
young systems should then fall somewhere in between the 
two histograms shown on Fig. ^1 Still, typical observed 
velocities will rarely exceed the order of 200-250 km s~^. 
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Fig. 10. Velocity distribution of TDGs with respect to the 
mass center of the merging pair of massive galaxies: real 
(3-D) velocity (green dashed) and projected along the line- 
of-sight, in absolute value (blue solid). Data accumulated 
from t = 500 to 2000 Myr, for TDGs surviving at least 
1 Gyr. 



3.3.4. Dark matter content 

We have measured the amount of dark matter inside the 
visible radius^ for tidal objects that survive more than 
1 Gyr. We find that inside this radius, the dark-to-visible 
mass ratio is on the average 0.09. It is never more than 
0.15, with the highest dark matter fractions generally 
found in the most massive TDGs. For comparison, massive 
spiral progenitors in our initial conditions have a dark-to- 
visible mass ratio of 0.65 inside their visible radius. 

Hence, tidal dwarf galaxies are strongly defi- 
cient in dark mat t er, ac cording to our simulations. 
iBarnes fc HernauTs^ l|l992|) found a similar result; our 
tidal dwarfs contain slightly more dark matter because 
we use more extended dark halos around the massive pro- 
genitors, so that more dark matter is found in the vicinity 
of forming TDGs. 



4. Discussion: comparison to observed TDGs and 
satellite galaxies 

4.1. Identifying young forming Tidal Dwarf Galaxies in 
colliding galaxies 

A large variety of star-forming entities have been re- 
ported in the tidal tails of real interacting systems: 
Super Star Cluster (SSC s), Gian t HII Complexes 
(GHC s 'l llGallagher et all [200lt iKnierman et alJ 



pnoi IWeilbacher et a] 



n 



2nn3albL ISaviane et a,].l l2004 

04: ilglesias-Paramo fc Vflchej 
12001b ICortese et al.ll2006j) and the more massive so-called 
"Tidal Dwarf Galaxies candidates" which, in fact, also 
correspond to a diversity of objects. However, several 
observed interacting systems are distinguished by the 
massive accumulations of matter, up to a few 10^ -^0' 
found at the tip of their tidal tails (e.g. IDuc fc Mirabel 
11994 Il998l: iHibbard fc Earned 12004 . In principle, they 
could be fake condensations caused by the projection of 
material along the line-of-sight in tidal tai ls seen edge- 
on. H owever their molecular gas content l(Braine et alJ 
and their large-scale kinematics provide evidences 
that they are real massive objects, while their internal 
kinematics, when resolved, i ndicates that these str uctures 
are gravitationally bound ijBournaud et al.l l2004l) . Most 
of these tidal objects are observed at a time when they 
are probably still forming, a few hundred Myr after the 
beginning of the collision. Without any detailed observa- 
tions and the help of evolutionary models, it is difficult to 
predict how they will evolve, and in particular whether 
they will become independent, stable, objects that would 
deserve to be called Tidal Dwarf Galaxies. Our large set 
of simulations show that among the objects formed at 
the base and the middle of the tidal tails, many are on 
rather high eccentric orbits which hampers considerably 
their ability to survive. Actually only a small fraction of 
them are still present after 500 Myr. About 75 percent 



^ defined here as the radius containing 90% of the baryonic 
mass (stars-|-gas) of the object 
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of the objects more massive than 10^ Mq - our mass 
threshold - survive less than 800 Myr. Such objects may 
be considered as TDG-candidates in on-going mergers, 
but will most likely be short-lived and not become dwarf 
galaxy satellites orbiting around the merger rem nant for 
at least a few dynamical times. iKroup studied 
the evolution of purely stellar low-mass satellite galaxies 
without dark matter (as expected for TDGs) in the 
potential well of the Galactic dark halo. He found that 
satellites on high eccentric orbits lose 10-20 % of their 
mass during each perigalactic passage. They nevertheless 
manage to survive reaching a stationary regime but their 
remnants contain then only 1 % of their in i tial m ass, 
about 10^ Mq in the simulations of iKroupal l)l997j) . To 
have final masses typical of dwarf spheroidals - the less 
massive known galaxies -, they then should initially have 
had stellar masses as high as 10^ — 10^" Mq, which can 
correspond only to the most massive objects, formed at 
the tip of tidal tails in our simulations. 

According to our simulations, the only objects that ac- 
tually become genuine dwarf galaxies are those formed at 
large radii, and preferentially near the extremity of the 
tails (Fig. [5J) via the mechanism studied by D04 which 
works for extended dark matter haloes. We then conclude 
that among the many structures observed in the tidal tails 
of young interacting systems, only the most massive ones 
that are formed near the tip of tails are likely to be the 
progenitors of Tidal Dwar f Galaxies, following the d efini- 
tion given e.g. by D04 or iHibbard fc Earned l)2004j) that 
requires them to be long-lived. 

One should however keep in mind the restrictions of 
these simulations. The produced TDGs were found to be 
long-lived essentially because they were orbiting on safe 
orbits (at large radii and low-eccentricities) and were not 
affected by a strong tidal field that could otherwise disrupt 
them. However we have not considered here the internal 
processes that could cause severe disturbances. Intensive 
feedback following an internal starburst may remove much 
gas from the TDGs ( Dekel fc Woo 2003 ). Massive objects 
may also fragment into smaller ones in higher resolution 
models. We will present in a separate paper high resolu- 
tion simulations zooming in onto TDGs that allowed us 
to conclude that, even taking into account their internal 
evolution, they are still long-lived objects that conserve 
most of their initial mass. 



4.2. Identifying Tidal Dwarf Galaxies around merger 
remnants 

Young forming Tidal Dwarf Galaxies are rather easily pin- 
pointed as long as they are still linked to their parent 
galaxies by a tidal tail. As explained earlier, an HI accu- 
mulation or CO detection at its tip make a it a promising 
TDG candidate. However, once they become "isolated" 
and that the umbilical cord is broken - typically tidal tails 
dissolve in 300-500 Myr -, they become much more dif- 



ficult to identify. Indeed, they resemble satellite galaxies 
orbiting the merger remnant. 

Several criteria to establish a tidal origin were dis- 
cussed in the literature: an unusually high metallic- 
ity received in in-he ritage from their parent galaxies 
iIDuc & Mirabe]lll998l). a deficiencv of non-baryonic dark 
matter (jBarnes fc Hernauistlll99 and thus a particular 
location in the TuUy-Fischer diagram (Hunter et al. 20o3). 
Several quests for identifying old TDG s among the dwar f 
galaxy population were carried out Csee lDuc et alJ l)2004bl) 
and references therein) , but sofar no dwarf satellite formed 
more than 1 Gyr ago during a collision has ever been un- 
ambiguously recognized, partly because the above men- 
tioned criteria are often difficult to check. 

The simulations presented here provide further and 
somehow simpler constraints to probe a tidal origin: 

Mass criterion The mass of TDGs is typically 0.1 to 1% 
the baryonic (gasH-star) mass of their progenitor when the 
latter had a gas content typical of bright spiral galaxies 
before the merger. It can be at the very most a few per- 
cent in the case of gas-rich progenitors. Even more mas- 
sive tidal galaxies were only able to form at high redshift 
when the parent spiral galaxies contained more gas. It is 
likely that such objects have already disappeared at red- 
shift 0. They would be anyway difficult to identify from 
their metallicity which is not expected to be very high. 
We can apply the mass criterion to some nearby systems 
(though still young) where the presence of massive TDGs 
had been suspected but later debated. For instance, the 
northern TDG-candidate of Arp 105 (Due et al. 1997) has 
a mass of nearly 10 percent of the parent galaxy mass 
(7 X 10^ solar masses in HI, the dominant component). 
An object with such a mass is then unlikely to be a real 
tidal dwarf, unless its mass had been over-estimated due, 
for instance, to the contribution of projected foreground 
or background material along the line-of-sight. The pres- 
ence of such a projection effect in the northern tail of 
Arp 105 is indeed suggested by th e large-sale kinematics 
of the tail (|^urnaud et aLll20n4|) The Southern TDG- 
candidate identified bv lDuc et al. l)l997|) . with a mass of 
0.5% of its parent spiral, is much more likely to be a gen- 
uine tidal dwarf. 

This mass criterion also suggests that the potential 
TDG-candidate NGC 3077 i n the g roup of M 81, linked 
to M 81 by a gaseous bridge l|Yun et al. 1994) is unlikely 
to be of tidal origin: according to its infrared luminosity 
and estimated HI mass, the mass of this galaxy is indeed 
as high as 15-20% of the mass of M 81. 

Spatial distribution Statistically, the simulated TDGs tend 
to be concentrated towards the equatorial plane of their 
massive host galaxy. However, this property is not a strong 
one for individual cases, since one given TDG can be found 
in a polar plane (see Fig.|SJ). Moreover, in the case of equal- 
mass mergers, the equatorial plane of the spiral progenitor 
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may hardly be observationally determined. Additionally, 
TDGs are rarely found at more than 15 times the op- 
tical radius of their progenitor, i.e. typically 200-250 kpc 
around a bright early-type host galaxy. In compact galaxy 
groups, the interaction with other galaxies may kick them 
from the host potential, and TDGs may be found further 
out from their progenitor. 

Velocity The relative velocities of TDGs, with respect to 
their hosts, are of the order of the circular velocities in 
dark haloes. Projected along the line-of-sight, this rarely 
results in velocities larger than 200-250 km s^^. Higher- 
velocity dwarf satellites are unlikely to be of tidal origin. 

Dark matter content Dwarf galaxies of tidal origin are ex- 
pected to be deficient in d ark matter (see Sect. 13.3."^ and 
iBarnes fc HernaiiistI 1 19921 )). Their dark matter content is 
negligible compared to their progenitor and to the many 
dwarfs with measured high M/L ratios. Note however that 
this prediction is made within a pure CDM frame. The 
presence of a baryonic and dissipative d ark matter com- 
ponent (e.g. IPfenniger fc GombesI ll99J) could radically 
change the result, for this component could participate 
to the formation of TDGs. 

Properties of the host galaxy Long-lived TDGs are formed 
during galaxy collisions in the range of mass ratio 1:1 to 
8:1, and more effici ently between 1:1 and 4:1. According to 
recen t simulations I Bournaud et alJ2005HNaab fc Truiillol 
I2OO5I) ■ the remnant of such barely-equal mass merger 
should have the morphology of an Sa or more likely an SO 
or even that of an elliptical, unless the colliding galaxies 
were particularly gas rich, as in the Early Universe. Tidal 
satellites should then be found more frequently around 
early- type galaxies. 

We did not explore parameters for which the parent 
galaxies do not merge and hence may remain spirals. In 
such fly-bys, which generally correspond to large impact 
parameters, the corotation radius is large; less material 
gets into tidal tails and the latter are shorter. If TDGs 
are produced at all, they will have lower masses. Indeed, 
the high- velocity or large-distance cases that we have sim- 
ulated (that are the closest to fly-bys without mergers), 
correspond to less numerous and less massive tidal dwarfs 
than the average (see Appendix A). 



4.3. Satellite galaxies around their hosts 

We have shown in the Section 3 that the distribution 
of TDGs around their progenitor galaxies is anisotropic. 
However, the projection on the sky may hide this property. 
Indeed, a TDG that in the 3D space lies in the equatorial 
plane of its massive progenitor, can be observed in the 
2D projection aligned with the minor axis of its massive 
host, causing possible confusion with a dwarf satellite on 
a polar orbit. 



We have statistically estimated the apparent angle 
between the host galaxy major axis and the apparent po- 
sition of the TDG, defined on Fig. This is the usual 
observational method used to study the dis tribution of 
dwar f satellites around bright galaxies (e.g. lYang et all 
I2OO6I and references therein). To infer the distribution of 
from that of the latitude i in three dimensions, we assumed 
that the host-l-TDG systems are observed under isotropi- 
cally distributed line-of-sights. The calculation was made 
analytically rather than assuming a limited number of pos- 
sible line-of-sights in simulations. The result is shown of 
Fig. nil : after projection on the sky plane, the anisotropy 
in the distribution of TDGs is still visible. 

Most studies of the distribution of dwarf satellites 
around massive host galaxies focus on the mos t mas- 
sive dwarfs. For instance, the work of iBrainerdI l|2005fl 
based on SDSS data includes only satellites more mas- 
sive than 1/8 of their host mass. Such objects are unlikely 
to be of tidal origin, since they are 1-2 orders of mag- 
nitude above the typical masses of long-lived TDGs. A 
st udy including lowe r mass objects has been performed 
bv lYang et alJ l(2006() . These authors found an anisotropic 
distribution of dwarf satellites, similar to what is found 
the simulations tidal dwarfs. These results cannot be di- 
rectly compared, since certainly not all the dwarfs are of 
tidal origin (see next section). However, TDGs may sig- 
nificantly contribute to forster the observed anisotropy in 
the distribution of dwarf satellites. Hence, the anisotropy 
of dwarf satellites may not only relate to the shape of 
dark haloes, but also to the tidal origin of some satel- 
lites. Therefore an anisotropic distribution could in prin- 
ciple be present eve n within spherical dark halos. Also 
note that lYang et alJ l|200d|) find that the distribution of 
satellites is more anisotropic around red, early-type host 
galaxies. This observational fact could be accounted for by 
the anisotropy of TDGs since tidal dwarfs are more likely 
to be found around early- type hosts. 



4.4. Contributions of TDGs to the dwarf population 

The rate of TDG formation in galaxy mergers cannot be 
directly deduced from our set of simulations, for we de- 
liberately privileged favorable cases: the distribution of 
mass ratios and orbital parameters in our simulations is 
not aimed at being cosmologically realistic. However, we 
have derived the constraints on parameters for the forma- 
tion of long-lived tidal dwarfs (Section 3.2) and this can 
be used to estimate the actual frequency at which TDGs 
are formed. 

Our simulations indicate that on average a bright spi- 
ral galaxy produces 1.3 tidal dwarf (considering only those 
that survive at least 1 Gyr), whenever it undergoes a 
merger with a galaxy between 1/4 and 8 times its mass, on 
a prograde orbit, with an orbital plane inclined up to 40 
degrees with respect to its disk plane. Assuming that real 
orbits for galactic encounters have an isotropic distribu- 
tion, this corresponds to a production rate of 0.4 TDG per 



10 



Bournaud & Due: From Tidal Dwarf Galaxies to Satellite Galaxies 



1.2 



Simulated TDGs (average 42.5°) 
SDSS (average 42.2°) 
Isotropic 




20 40 60 80 

Inclination (projected) 

Fig. 11. Distribution of the projected inclination of dwarf 
satellites of tidal origin w.r.t. their host (progenitors) 
galaxy plane: the angle 6 between the dwarf galaxy and the 
host galaxy apparent major axis is measured (see Fig. [71). 
The distribution has been derived assuming that the sys- 
tems are randomly projected on the sky. The average value 
of is 42.5 degrees, while an isotropic distribution corre- 
sponds to = 45 degrees. The shghtly more anisotropic 
distribution observed in the Sloan Digitized Sky Survey 
by Yang et al. (2006) is also shown. 



merging spiral. For mergers in the 1:1-4:1 range of mass 
ratios, both spirals can form long-lived TDGs, and the 
production rate increases to 0.8 TDG per merging pair. 

This estimate a priori assumes that the distribution of 
impact parameters and initial velocities in our sample of 
simulation is statistically representative of real colliding 
systems, which is probably not the case. However, our 
results summarized in Table 1 show these two parameters 
have a minor influence on the formation of TDGs. Then, 
even if their distribution in our sample is not realistic, this 
does not result in an important bias on the estimated rate 
of TDG production. 

lOkazaki fc Tanieuchil 1)2000 *) claimed that all dwarfs 
in the Universe could have a tidal origin. Their estimate 
relied however on a production rate of 1 to 2 tidal dwarfs 
having a live time of at least 10 Gyr, which is actually 
much higher than that predicted by our simulations. 

The number of TDGs versus time in our simulations 
(207 at t = 500 Myr, 143 at i = 1 Gyr, and 119 at i = 
2 Gyr) fits well with an exponential decay and a lifetime 
of 2.5 Gyr, so that only 20 percent of the TDGs actually 
survive up to 10 Gyr. Therefore, our best production rate - 
0.8 TDG with life of 1 Gyr per merging pair - corresponds 
to just 0.1-0.2 TDG with a life of 10 Gyr, ten times lowe r 
than the value assumed bv lOkazaki fc Taniguchil (|2D0rf). 
Using their other hypotheses, in particular the evolution 
of the number of mergers with redshift, we conclude that 
at most 10% of the dwarf galaxy population could be of 
tidal origin. 

In fact, the study bv lOkazaki fc Taniguchil l)2000|) im- 
plicitly covers a larger range of mass ratios. Only about 
one third of the mergers are likely to fall within the 1:1- 



8:1 range favorable to long-lived TDGs. Since the latter 
estimate is very tentative, we will simply conclude that a 
few percent of the dwarf galaxy population can be of tidal 
origin. On the other hand, we did not take into account 
the fact that more distant galaxies were more gas rich and 
had probably an HI disk much more extended than in lo- 
cal spirals, two properties which may enhance the TDG 
production rate. 

Only few observational studies have sofar tried 
to quantify the contribution of old, "detached" 
(thus long-lived) TDGs to the dwarf population. 
iDelgado-Donate et alJ l)2003j) did not found any excess of 
dwarfs around a sample of strongly interacting systems 
and concluded that the contribution of long-lived TDGs 
should be rather minor. On the other hand, in denser 
environments, multiple collisions could kick tidal dwarfs 
from the progenitor potential and thus contribute to 
increase their life expectancy. iH^msberger et al.l ljl99fj) 
probed the dwarf population of a sample of compact 
groups of galaxies and concluded that 50% of them 
were probably tidal objects. The production of tidal 
dwarfs in the cluster environment could be hampered 
by the fact that tidal tails develop in a less efficient 
way a nd tend to b e more diffuse in massive potential 
wells l|Mihosl l2004fl . However such objects may have 
formed in sub-group s falling in the cluster. A rp 10 5 
( Due fc M irabel 1994'!. NGC 5291 llDuc fc Mirabe]lll993) , 
IC 1182 (,Iglesias-Paramo et all l2003l) are instances of 
systems belonging to clusters or sub-structures within 
clusters where young TDG candidates have actually been 
identified. The presence of a large popu lation of cluster 
dwar f ellipticals in Com a and Fornax I Poggianti et alJ 
l200lt iRakos et al.1l2000l). of dwa rf irregulars in Hydra-I 
and Hercules ||Diic et all l200ll llglesias-Paramo et alJ 
1^03) having deviant - high - metallicities for their 
luminosity could be interpreted as hints for a tidal origin 
for some of them alt hough other mechan isms could 
explain these anomalies i Conselice et al.lf2003|) . 



4.5. Looking for TDGs in the Local Group 

The distribution of the dwarf satellites of the Milky- 
Way is obser ved to be anisotropic within so-called 
Great Plane(s) l)Lvnden- BcU 1976; Kunkcl 1979; Mai ewskij 
19941). This fact was expl oited bv lKrouoa et al.l |2005h and 
Sawa fc Fuiimotol ^^2005^ to claim that the dSphs around 



the MW could be of tidal origin. However, the satellites 
are actually not concentrated towards the disk plane, but 
along one or two nearly perpendicular directions. Our sim- 
ulations indicate that such a distribution is a priori in 
confiict with the hypothesis that the M-W dSphs are old 
TDGs since the latter are statistically concentrated to- 
wards the equatorial plane of their progenitor^ ... unless 
the parent galaxy is not the M-W. Let us imagine that long 
ago a spiral with a disk perpendicular to the M-W disk, or- 
bited in a polar plane w.r.t the M-W and produced in that 



though high latitudes TDGs can also be present; see Fig|Sl 
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plane long-lived TDGs. Nowadays, we do not see the disk 
of this progenitor, either because it has flown far away 
or it has merged with the M-W disk. According to our 
simulations, the large number of putative TDGs - more 
than 10 dSphs are orbiting around the M-W - indicate 
that either the collision at their origin occurred less than 
1 Gyr ago, before most of them are destroyed - however, 
no sign of such recent encounter is clearly visible, except 
the LMC/SMC/interaction -, or that several mergers, all 
with specific orbits included in the Great Plane(s), oc- 
curred long ago and produced each 1-2 long-lived TDGs. 
However, the present morphology of the Galaxy is incon- 
sistent with the idea that it has merged in the past with 
several companions with mass ratios of at least 8:1, the 
limits we found to produce a TDG: if this had occurred, 
the M-W would hav e become an SO or E galaxy (e.g. 
iBournaud et al.ll2005|) . 

In that respect, M 31 is more likely to host tidal 
satellites; its morphology has an earlier type and nu- 
merous signs of possible past interactions, such as long 
stellar strea i ns, w er e found in its outer most regions 
l|lbata et all l200l|) . iKoch fc Grebe] l|2006^ discuss the 
three-dimensional location of the satellite galaxies around 
M 31. The global distribution of satellites in this system 
do not show any clear concentration towards a particular 
plane, but the early-type dwarfs and the dSphs in particu- 
lar seem to be concentrated towards a polar plane. As for 
the case of the Milky- Way, this distribution would require 
that several galaxy mergers have occurred on this config- 
uration, which makes the tidal scenario for the origin of 
M 31 's dwarfs unlikely. 

Therefore, just from the distribution and number of 
the dSphs located in the vicinity of the M-W, and the sta- 
tistical results obtained with our simulations, one should 
conclude that it is unlikely that the majority of them are 
of tidal origin. This is also the case for M 31. Other sce- 
narios to expla i n the i r distribution in g reat planes (e.g. 
iKoch fc Grebeil l2006t iLee fc Kan3 12006(1 are then more 
realistic. 

However, other, somehow more direct, criteria to probe 
a tidal origin should be investigated. For insta nce, a devi- 
ation from the luminosity-metallicity relation ( Due et al\ 
I2OOQ1') due to prc-cnrichment processes in the parent's 
gaseous disk or the absence of a non baryonic dark mat- 
ter component (see Sect. 13.3^ are clear characteristics of 
TDGs. 

In that respect, the Local Group dwarfs were subject 
of active debates. They are usually considered as being 
the objects showing the highest M/L ratios (up to 100) 
and thus containing the highest quantities of dark matter, 
which would argue against a tidal origin, at least in a CDM 
context. Nevertheless the accurateness of the estimate of 
their total mass is controversial, and claims were made 
that it had been largely over-estimated due to the lack of 
statistics - the dynamical mass is derived from the veloc- 
ity dispersion of a limited number of individual stars - o r 
due to the tidal shaking of the stars ijKroupa et alJ Eoo3) . 
Beside, the Local Group dwarfs are remarkable for their 



variety of M/L ratios, chemical enrichment and more glob- 
ally star for mation history despite being all in the same 
environment l|Grebe]ll200il . This leaves room for a tidal 
origin in some of them. 

4. 6. A note on merger induced globular clusters 

In our analysis, we only considered sub-structures with 
masses exceeding 10^ Mq to make sure they contained 
enough particles and had a fair resolution. Our numerical 
study did not aim at making predictions on lower mass ob- 
jects, in particular the potentially forming globular clus- 
ters^. Objects of masses typically 10^ — 10^ Mq could be 
numerous in tidal tails, and their survival less problematic 
than the objects with typical masses of 10* Mq. Indeed, if 
concentrated enough, they will be less affected by tidal dis- 
ruption, unde rgo less dynarnical f ri ction and face a slower 
orbital decay l|Miocchi et alJ2006|) . lKrour)all)l997|) studied 
the evolution of 10^ Mq dark matter free satellites and 
concluded that they could be long-lived, as well as globu- 
lar clusters possibly formed in mergers. Higher resolution 
simulations including a description of the formation pro- 
cess are needed to study this scenario of globular cluster 
formation in mergers. 



5. Conclusion 

In this paper, we have studied the formation and the evo- 
lution of Tidal Dwarf Galaxies and related candidates, 
in about 100 numerical simulations of collisions of spi- 
ral galaxies, with gas content typical of gas-rich low to 
moderate redshift spirals. The identification of massive 
{M > 10* Mq) substructures formed in tidal tails, and 
their follow-up during 2 Gyr, have shown that: 

— Among the numerous tidal objects, those suscepti- 
ble to survive more than 1 Gyr were formed in the 
outer parts of tidal tails; they also turn out to be the 
most massive ones. Gravitational clumps which grow 
all along the tails rarely survive more than 500 Myr, 
mainly because they are formed on more eccentric or- 
bits, and generally fall back onto their progenitors or 
are disrupted by the tidal field losing much of their ini- 
tial mass. These objects should not be considered as 
genui ne Tidal Dwarf Gal axi es, according to the defini- 
tion oflDuc et alJ ^2004b|l orlHibbard fc Rarnej HqoI , 
even if they have, at the moment they are observed, 
the apparent characteristics of young forming dwarf 
galaxies. 

— Long-lived tidal dwarf galaxies are preferentially 
formed on prograde, close-to-coplanar orbits (inclina- 
tion smaller than ~ 40 degrees). In such cases, more 
than 2 long-lived TDGs are rarely formed from one 
progenitor galaxy. These results have led us to the con- 
clusion that a small but significant fraction of a few 



Note though on Fig. 1 the presence of low mass objects in 
the vicinity of the merger remnant. 
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percent of dwarf satellites could be of tidal origin. The 
fraction of tidal dwarfs can be larger around early-type 
hosts that have undergone numerous mergers or in rich 
environments where tidal dwarfs can be longer-lived. 

— The statistical properties of the long-lived TDGs, 
those that become dwarf satellites, have been studied. 
They are typically found within 10-15 optical radii 
from their progenitor, with moderate relative veloc- 
ities, and masses at the very most a few percent of 
their progenitor mass (or about 10% of their progen- 
itor initial gaseous mass). Their spatial distribution 
is anisotropic, flattened towards the equatorial plane 
of the parent spiral galaxy. Would the contribution of 
TDGs be important in some environments, they would 
then contribute to foster the anisotropy already ob- 
served for the distribution of the general dwarf satellite 
population. 

— Given their number and spatial distribution in po- 
lar Great Planes, the majority of the satellites of the 
Milky- Way and M 31 are unlikely to all be of tidal 
origin. 

These results can help identifying, in real young inter- 
acting systems, the tidal objects that are the most likely 
to be long-lived and provide hints on the location of the 
long searched "old" TDGs which, once the tidal tails link- 
ing them to their parents have disappeared, are difhcult to 
identify. If searched randomly without any selection crite- 
ria, tens of dwarf satellites would have to be analyzed for 
a few TDGs to be found. 

Galaxy mergers can also produce a significant number 
of smaller mass objects like globular clusters, but the reso- 
lution of the present simulations did not enable us to study 
them in detail. Also, the internal properties of TDGs have 
not been studied in this paper. They will be the subject of 
a forthcoming study using higher resolution simulations, 
dedicated to resolve the internal dynamics of tidal sub- 
structures, in particular those that are now known to be 
long-lived. 
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Appendix A: Parameters and results of numerical 
runs 

We give in Table lA.ll the parameters used for each sim- 
ulation. The mass ratio is M . The orientation of the or- 
bit, prograde (P) or retrograde (R) is given for the most 
massive galaxy and the less massive one respectively. The 
impact parameter R and relative velocity at an infinite 
distance V are computed assuming that dynamical fric- 
tion is negligible before the beginning of the simulation. 
The initial inclination of the orbital plane to the disk plane 
is given by ii and 12 for each colliding galaxy. iVi and N2 
are the number of long-lived tidal dwarf galaxies (surviv- 
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Table A.l. Run parameters and results: number of long-lived TDGs (surviving at least up to f = 1 Gyr) formed in 
the material of each parent galaxy. Parameters are defined in text. 
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Table A. 2. Run parameters (continued). 
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Fig. 1. Face-on snapshots for run 69 (coplanar encounter). Red and blue colors correspond to regions dominated by 
old stars vs. gas and young stars. The system is shown at four different instants : (a) at t = 150 Myr after the 
pericenter, gas accumulation progenitors of TDGs A and B are already visible at the tip of the tails - (b) Location of 
the 7 TDG-candidates more massive than 10^ Mq. Initial masses (in 10® Mq) A: 14, B: 8, C: 1.5, D: 2.4, E: 7, F: 6, 
G: 1.2. - (c) Sitt — 400, object G has fallen back onto the progenitors, and objects D, C, E, F begin to fall-back, while 
objects A and B remain at larger radii. - (d) at i = 800, only objects A and B still exist according to our criterion. 
Remnants of structures F and D are visible, but they have lost much of their mass: they now have masses of 0.6 and 
0.8 respectively (object F has passed within 20 kpc from the merger remnant center). TDGs A and B, located at larger 
radii, did not undergo such a mass loss: they now have masses of 12 and 6 respectively. Object B falls back onto the 
merger remnant at f = 1.7 Gyr; object A survives up to the end of the simulation at i = 2 Gyr. 



Bournaud & Due: From Tidal Dwarf Galaxies to Satellite Galaxies 



17 



90 



05 60 



30 



TDGs 

Isotropic 

R < 82 kpc 



90 



05 60 



30 



TDGs 

Isotropic 

82 kpc < R < 121 kpc 



20 



40 60 
Inclination 



80 



20 



40 60 
Inclination 



80 



90 



m 60 
j3 



TDGs 

Isotropic 

121kpc<R< 157 kpc 



30 



90 



05 60 

J2 



30 



TDGs 

Isotropic 

157 kpc< R 



I < \ < \ < \ < I I , \ , \ , \ I \ I 

20 40 60 80 20 40 60 80 

Inclination Inclination 

Fig. 6. Same as Fig. but for long-lived TDGs separated into four radial bins of equal population (radial quartiles 
at 28, 121, and 157 kpc). The anisotropy of the distribution of long-lived TDGs increases with radius. 




Fig. 7. Definition of the real inclination of a TDG with respect to the disk plane (angle i in the 3-D space, left), and 
its apparent position w-r-t the projected major axis of its host (angle 9 in the 2-D projection plane, right). One given 
object at one given instant has a single latitude i but a probability distribution for 6 varying from to 90 degrees, 
depending on the line-of-sight. 



